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Motivated by the quasi-lD antiferromagnet CaV204, we explore spin-orbital systems in which the 
spin modes are gapped but orbitals are near a macroscopically degenerate classical transition. Within 
a simplified model we show that gapless orbital liquid phases possessing power-law correlations may 
occur without the strict condition of a continuous orbital symmetry. For the model proposed for 
CaV2 04, we find that an orbital phase with coexisting order parameters emerges from a multicritical 
point. The effective orbital model consists of zigzag-coupled transverse field Ising chains. The 
corresponding global phase diagram is constructed using field theory methods and analyzed near 
the multicritical point with the aid of an exact solution of a zigzag XXZ model. 
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I. INTRODUCTION 

The orbital degree of freedom plays a major role in 
the physical properties of strongly correlated transition 
metal oxides.^ For Mott insulators in which the low en- 
ergy electrons occupy orbitals with Cg symmetry, or- 
bital order typically occurs due to Jahn- Teller distor- 
tions at relatively high temperatures and determines the 
spatial anisotropy of spin exchange interactions at low 
temperatures.l^l However, in t2g systems the Jahn- Teller 
interactions are much weaker and the orbitals can remain 
unquenched down to low temperatures.'^ The stabiliza- 
tion of orbital phases may eventually depend on the cou- 
pling between spins and orbitals. This coupling can stem 
either from the on-site relativistic term or from Kugel- 
Khomskii^- type coupling. Particularly intriguing is the 
possibility of quantum melting of orbital order. One- 
dimensional (ID) models with SU(2) symmetry for both 
S — 1/2 spins and t — 1/2 orbital pseudospins are known 
to host liquid phases.!^ On the other hand, the ID model 
with S = 1 spins accounts for the orbital Peierls state 
observed in YVOa.'^In three dimensions an orbital li(^ 
uid (OL) was proposed as the ground state of LaTiOa,^ 
but was later dismissed by experiments.13 

In fact, it may be even possible to invert the usual or- 
der and freeze out the spins at higher energies than the 
orbitals. Here we will show that this can be achieved 
in ID Haldane gapped spin systems in the vicinity of a 
classical orbital transition. Our theory is motivated by 
experiments on the S = 1 zigzag compound CaV204.'^ 
In this material the twofold single-ion orbital degener- 
acy is lifted by a structural transition below Ts « 141 
K, whereas the strongest intrachain exchange coupling is 
estimated around J « 230 K. Antiferromagnetic order 
only sets in at Tn ~ 71 K. Experiments favor a scenario 
of orbital ordering that turns the S = 1 zigzag chains 
into S ~ 1 ladders .'Si Recently, Chern and Perkind^ 
introduced a model for CaV204 which contains Ising 
antiferro-orbital (AFO) interactions. These interactions 
stabilize a classical AFO ordered phase at weak struc- 
tural distortion (small orbital energy splitting) , whereas a 



para-orbital (PO) phase corresponding to the experimen- 
tally favored ladder structure occurs at strong distortion. 
Their model at zero structural distortion already presents 
rich physics occurring as function of increasing relativis- 
tic on-site spin-orbit coupling: Chern et aZ.I^ found that 
first an Ising transition turns the spin 5 = 1 system from 
the Haldane spin-liquid to a Neel phase, across which the 
orbitals maintain their finite Ising order, which however 
is destroyed at a second Ising transition at yet stronger 
spin-orbit coupling (they have considered ferro-orbital 
rather than AFO interactions, which however leads to 
similar physics). More recently, Nersesyan et al.^^ found 
that whereas the above two-stage ordering scenario oc- 
curs for the limiting case of large orbital gap as compared 
to the spin gap, in the opposite limiting case of large 
spin gap a single Gaussian transition between the Hal- 
dane spin liquid and the Neel spin-orbital ordered phase 
takes place as a function of the strength of the spin-orbit 
coupling. In close connection, spin-orbit interaction ap- 
pears to be important for the still unclear orbital order 
in related spinel compounds.l^^lliS 

In this work we consider the effects of weak spin-orbit 
interactions near a macroscopically degenerate classical 
orbital transition, such as the AFO to PO transition 
which occurs in the model proposed for CaV204.'^ We 
show that this degeneracy renders the orbitals suscepti- 
ble to arbitrarily small quantum fluctuations induced by 
coupling to gapped spin modes. As a result, the clas- 
sical transition point develops into a multicritical •point 
(MCP). We describe the phases that appear near the 
MCP and the universality classes of the transitions. 

The outline of the paper is as follows. In Sec. [II] we 
study a simplified model in which an Ising orbital chain 
and a Haldane spin chain are coupled only by relativis- 
tic spin-orbit interaction. We show that by integrat- 
ing out the spins one obtains an effective orbital model 
whose phase diagram contains a MCP, with two inter- 
mediate phases induced by the spin-orbit coupling. One 
of these phases is a critical orbital liquid, described as 
a Luttinger liquid (LL) with power-law correlations, and 
forms without the strict condition of continuous orbital 
symmetry. In Sec. [III[ we study the model proposed to 
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describe CaV204, which besides relativistic spin-orbit 
coupling contains Kugel-Khomskii type interactions. We 
consider particularly the experimentally relevant regime 
in which the next-nearest-neighbor exchange coupling in 
the zigzag structure is much stronger than the nearest- 
neighbour coupling. In order to construct the phase di- 
agram of this zigzag model, we first present an exact 
solution of a U(1)®U(1) symmetric model of hard core 
particles in Sec. |IV[ We find that the phase diagram of 
the latter also contains a MCP, but with a single inter- 
mediate gapless phase interpreted as an orbital Luttinger 
liquid that breaks the symmetry between the legs of the 
zigzag chain. We then return to the model for CaV204 
and show that in this case the orbital Luttinger liquid is 
unstable and gets replaced by a gapped phase with co- 
existing order parameters. The existence of this phase in 
the limit of weak spin-orbit coupling is supported by a 
global analysis of the phase diagram based on field theory 
methods carried out in Sec. |Vj Finally, the conclusions 
are presented in Sec. |VI[ 

II. TOY MODEL 

We start by considering a spin-orbital toy model con- 
taining on-site relativistic spin-orbit coupling: 

H = YS^i^Ui - K + J^J ■ S.+i - SArJ^I). (1) 

Here t is an orbital pseudospin-1/2 operator representing 
nearly degenerate t2g real orbitals (|t^ = i) — \yz) and 
= —\) = \xz^^ and S is a spin-1 operator. The 
spins interact via exchange coupling J and the orbitals 
via the Ising interaction A, which reflects the directional 
nature of the orbitals. Physically, both A and J originate 
from kinetic exchange processes in Mott insulators and 
are expected to be of the same order. The orbital field 
h depends on the strength of the structural distortion 
that lifts the orbital degeneracy. We consider /i as a 
parameter, controlled e.g. by pressure. The last term 
represents the projection of the spin-orbit interaction AL- 
S, where L is the orbital angular momentum, in the two- 
dimensional orbital subspace. 

A ferro-orbital version of model ([T]) was studied in 
Ref. jTI] Here we consider antiferro-orbital coupling 
A > 0, which arises naturally in the vanadates. liSIill xhe 
antiferro-orbital model with h — Q was studied in Ref. [T2| 
using field theory methods. For h — Q, it was found that 
in the regime of J 3> A, a Gaussian transition separates 
a Haldane spin liquid AFO phasd^ at small A from a Neel 
spin-orbital phase at large A. In the opposite regime of 
A ^ J, the Gaussian transition splits up into two Ising 
transitions where upon increasing A Neel spin order de- 
velops before the Neel order of orbitals disappears. In 
the following we study the model at finite h, which is 
interesting because the orbital energy splitting competes 
with the Ising interaction A. We also focus on the limit 



of weak spin-orbit interaction, in which the excitations 
can be regarded as predominantly orbital or spin, as op- 
posed to hybrid spin-orbital excitations. For A = 0, the 
spins are in the Haldane phase, which has an energy gap 
of about As « 0.41J,^^ and the orbitals are described as 
a classical Ising chain. The orbital sector exhibits a clas- 
sical phase transition at \h\ = A: for \h\ < A the ground 
state has Neel order of (AFOz state) and for \h\ > A 
it is a PO state with orbitals polarized in the direction 
of the field (hereafter we assume h> 0). The crucial as- 
pect of this transition is its macroscopic degeneracy: at 
the point h — A, X = 0, all states that do not contain 
nearest-neighbor (nn) down pseudospins are degenerate. 
The number of configurations satisfying this constraint 
scales like ip^ , where ip = (1 + \/5)/2 is the golden ratio 
and L is the system size.liSl 

The macroscopic degeneracy of the classical critical 
point should be lifted by arbitrarily weak quantum fiuc- 
tuations. In the regime A A, J the spin chain remains 
in the Haldane phase with a finite spin gap. On the 
other hand, for ft, « A the orbital gap becomes small. 
At energy scales far below the spin gap the spin chain is 
in a singlet ground state and the effective orbital model 
can be obtained by integrating out the spin fluctuations 
within second-order perturbation theory. 

The procedure of integrating out the spin modes in the 
limit A <C J was performed in Ref. J^in the absence of 
the orbital field. Here we do the same for finite h. We are 
mainly interested in the vicinity to the critical point, i.e. 
A~J^|A — /i|,A. Second order perturbation theory 
in the spin-orbit coupling in Eq. ([T]) yields an effective 
Ising orbital interaction-^l 

j e>i 

where ^ a"^^^^'' proportional to the spin-spin cor- 
relation in the Haldane ground state and decays expo- 
nentially for distances larger than the correlation length 
^ ~ 6. For a qualitative understanding of the resulting 
orbital phases it suffices to consider the nn term £ = 1 
as the leading transverse orbital coupling. The effective 
orbital model then becomes 

i/eff.. - E(^^/^+l - K + ^^/^jV^)' (3) 
j 

with £ = I and k ^ A^/As ^ A. We keep the integer 
index £ in Eq. (|3| in order to refer to a family of models 
which will be convenient later. We recognize -/?cff , i ^ 
the axial next-nearest-neighbor Ising (ANNNI) model in 
1 + 1 dimensions. It has been extensively studied in sta- 
tistical mechanics and finds applications, for instance, in 
the physics of submonolayer adsorbates. For a review see, 
e.g., Refs. I19I20I Based on various analytical and numer- 
ical methods it is known that the phase diagram of the 
ANNNI model contains a multicritical point (MCP) at 
h — A, K = 0, from which four phases emerge as shown 
in Fig. la. 
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Table I: Symmetry characterization of the orbital phases oc- 
curring in the ANNNI model Eq. ([sjl with 1=1. The ph ases 
are characterized in terms of the order parameters in Eq. Q. 
For an intuitive picture we include cartoons of the orbital 
configurations, with = |) = | t)i = — |) = I i)i 



pseudospin fluctuating in the xy plane. 



), and |») denotes a 



The phases are specified by broken symmetries with re- 
spect to Eq. (|3|, namely: (i) translation and (m) Ising Z2 
symmetry — r^, associated with order parameters 
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respectively. The PO phase at large h has a unique 
ground state and is disordered, 77^ = 77^ = 0, whereas 
the AFOz phase is doubly degenerate and is character- 
ized by finite M^; see Table |T] which also illustrates the 
pseudospin configuration. Two other phases appear be- 
tween AFOz and PO. The phase denoted AFOy, which 
is favored by large k, breaks the Z2 symmetry and has 
AFO order along y, hence finite (we assume k > 0). 
Generically, there cannot be a direct AFOy- AFOz tran- 
sition within a Ginzburg-Landau theory of two indepen- 
dent order parameters, unless extra symmetry is present. 
Indeed, as implied by this symmetry argument, a so 
called "floating phase"jhere interpreted as an orbital liq- 
uid (OL), is knowiPllsa to arise in between the two; see 
Fig. 1. In this gapless phase pseudospins fluctuate in 
the xy plane and have an emergent U(l) symmetry cor- 
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Figure 1: (a) Phase diagram of effective orbital model ([TJ near 
the multicritical point at /i = A, A = 0. Here k oc A^/A^. The 
properties of the four orbital phases are described in Table 
|l] (b) Schematic global phase diagram of the ANNNI model 
given in Eq. with i=\. 



responding to a conserved z-pseudospin magnetization. 

Note, however, that if we extrapolate model ([ij be- 
yond the perturbative regime, an explicit U(l) symme- 
try occurs along the line ft, = at k = A. At this point 
there is a direct transition from the AFOz to the AFOy 
phase. The global phase diagram of the ANNNI model 
is shown schematically in Fig. Ib.^^iHll! We believe that 
the U(l) symmetric point of the ANNNI model is in cor- 
respondence with the Gaussian transition discussed in 
Ref. 12 for /i = and A <C J, between the Haldane spin 
liquid AFOz phase at small A and the Neel spin-orbital 
phase at large A. Indeed, in our AFOy phase flnite Tf> 
orbital order implies a finite Neel order for the spins, 
(X]i(~l) ''S'f ) ^ 0, and the sign of this Neel order parame- 
ter is tied to the sign of rf such as to minimize the A term 
in Eq. ( ij , hence the AFOy phase corresponds to the Neel 
spin-orbital phase. An effective field theory was devised 
in order to access the various phases of the ANNNI model 
from the U(l) symmetric point k/A — 1 = /i/A — 0.1^ 
However, it remains unclear both from field theorji^Sl and 
from numerical studie^^H-v^rhether or not the orbital liquid 
phase occurs already for infinitesimal h from the Gaus- 
sian transition at k/A = 1. Its existence was however 
confirmed numerically for /i/A > 0.05.1^ On the other 
hand the region of the phase diagram near the MCP 
ft = A, K = 0, has a transparent description in terms 
of effective hardcore particles j^^ '^^'as will be discussed in 
the following subsection. 



A. Analysis of the ANNNI model close to the 
MCP 

In this section we briefly review field theory arguments 
that lead to the phase diagram of the ANNNI model in 
Fig. 1. Similar arguments will be used in the construc- 
tion of the phase diagram for the zigzag model in Sec. 

nil 

For the ANNNI model, it is possible to obtain the tran- 
sition lines near the MCP exactly. Following Rujan,'22Jwe 
consider a deformation of the model (with £ = 1), 



Hxxzj 



H 



xxz 



SH,, = 



) + 



^3 ^3+e> 



(5) 



The ANNNI model is obtained by setting e = 1. The 
model with e = is the XXZ model, which has a U(l) 
symmetry associated with the conservation of r^^j = 
J2j Tj ■ The XXZ model is exactly solvable by Bethe 
ansatz, but we shall follow a much simpler approach valid 
close to the MCP.^^ We will use the XXZ model as start- 
ing point and then analyze the effects of the e perturba- 
tion in Eq. (|5|. 

It is convenient to introduce Jordan- Wigner (JW) 
fermions by = ^ — fij, fij 
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c^(-l)^e'''^'='"'. Then the XXZ model becomes 



xxz.i = ^Ahjhj+i - fihj - ^(c]c^+i + h.c.) 



and the perturbation reads 



(6) 
(7) 



The chemical potential jj. ~ A — h vanishes at the MCP. 
Near this point the largest energy scale is A, which repre- 
sents a strong repulsion between fermions occupying nn 
sites. Interestingly, in the limit A/k — >■ oo, the XXZ 
model can be solved exactly without Bethe ansatz by 
mapping to free fermions with an exclusion constraintl^Sl 
The interaction is replaced by the hard-core constraint 
that no two fermions can occupy nn sites. Exactly at the 
MCP, all states which satisfy the nn constraint are de- 
generate and form a reduced low energy subspace. The 
pairing operator 5Hi vanishes within this subspace. To 
first order in k/A, the degeneracy is lifted by hopping of 
new fermions with annihilation operator aj on a reduced 
lattice with length Lf = L(l — v) where = /cf/tt is the 
filling fraction, described by the effective Hamiltonian^S- 



Hxxz,i = ["1 + "-l+i"-]) - ■ (^) 

3 = 1 

The ground state energy for fixed v is easily evaluated, 
egs{i^) = lini ^ 

Minimizing egs{i') with respect to i', one find^^H 



h K 
A^^^2A 



COSjt^ 1 . TTU 

— sm 
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(10) 



As a result, the MCP of the XXZ model (in the absence 
of SHi) is a triple point with a critical region in which 
1/ changes continuously in the phase diagram for — < 
^ — 1 < separating the AFOz with 1^ — 1/2 from 
the PO phase with v = 0. The low energy physics of this 
orbital liquid (OL) phase is described by the Luttinger 
modeP3 



HLL = \j dx [K{dj f + K-\d,<t>f] , 



(11) 



with [(l){x)^dx'0{x')\ = i5{x — x'), and the bosonization 



formula of JW fermions is c. 



c{x) 



Here v is the velocity, and K is the 
Luttinger parameter. The velocity v is determined by 
the energy spacing for a finite size system A^ = ^v. 



From Eq. (9 1 we find A^ 



.in (3^) 



givmg 



sm 



(iS)- Luttinger parameter is re- 



lated to the compressibility bjl^H K — ttv (J^t^ 

Using Eq. (l9|, one finds to first order in k/A that 
K = {1 — vyi^ The transition to the PO state in the 
XXZ model corresponds to the limit of vanishing den- 
sity, J/ — 0, where K ^ \. The transition to the AFOz 
state corresponds to half-filling v — )■ 1/2, in which case 
if — > 1/4. The latter is a commensurate-incommensurate 
(C-IC) transition.!^ 



The phase diagram of the ANNNI model can be in- 
terpreted in terms of the instability of the critical phase 
against (for general < e < 1). Within second 
order perturbation theory, besides fermion-fermion in- 
teractions, a next-nearest-neighbor (nnn) pairing term 
SHnnn ^ ('*^/^) X]j + h.c.) is generated in the 

reduced low energy subspace. In the bosonic language, 
the leading U(l) symmetry breaking perturbation is the 
operator cos(2-\/7r6').'^In our notation the latter has scal- 
ing dimension 1/K and becomes relevant when K > 1/2. 
Since ii' = 1/4 at the C-IC transition, 6 Hi is irrelevant 
near the boundary with the AFOz phase and the sys- 
tem has an emergent U(l) symmetry in the OL phase. 
Moving towards larger h, a Kosterlitz-Thouless transi- 
tion takes place at the line where K = 1/2, which can 
be found using Eq. mn at 1/ = 1 - I/V2. Beyond this 
line the U(l) symmetry breaking term 6H1 opens up a 
gap in the spectrum. The 9 field is pinned, corresponding 
to long-range order in the plane perpendicular to the z 
direction. For k > in Eq. this implies Neel order 
of Tj along the y direction, with order parameter ry^ 7^ 0. 
Finally, if we keep increasing h, an Ising transition into 
the polarized PO state takes place. 



We note that adding further neighbor interactions to 
Eq. (|3|, one can argue that as long as these can be 
treated as small perturbations, the phase diagram in Fig. 
[TJremains qualitatively correct. The perturbations renor- 
malize the parameters of the effective field theory and the 
transition lines can no longer be determined exactly, but 
the Luttinger parameter still assumes the universal value 
K — 1/423 at the transition to the AFOz phase, and 
the OL is stable against U(l) symmetry breaking per- 
turbations. However, when further neighbor interactions 
dominate, the phase diagram can change qualitatively. 
In any case, it is guaranteed by the macroscopic degen- 
eracy of the critical point that new phases will emerge 
above it. The case of dominant nnn interactions will be 
discussed in Sec. IIIII 
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III. ZIGZAG MODEL 



We now carry out the same program as in Sec. [ll|for 
the spin-orbital model for CaV20i^ 



HJiO, - J[0,)Sj ■ S,+i], (12) 

where Ji, J'l, J2, A > 0. Compared to the toy model 
in Eq. {[ij. here the nn spin exchange interaction de- 
pends on the orbital configuration (Kugel-Khomskii 
type interactions) via the antiferro and ferro-orbital 
bond operators Oj — 5(1 — '^TjTj^^) and Oj — 

+ (— 1)^tJ] [i + (— • In the zigzag structure 
the nnn interaction J2 arises due to direct exchange of 
occupied \xy) orbitals. Experiments^ on CaV204 suggest 
that J2 S> Ji ^ J(, implying that spins are close to the 
limit of weak coupling between even and odd sublattices. 

We start with symmetry considerations. Model (12) 



has the following discrete symmetries: (i) time reversal 
symmetry S — — S, r^'^ — > t^'^, — — r^; (ii) reflec- 
tion about either even or odd bond centers: j — > — j + 1 
or J — i' — j — 1. Note that for ft, > the nonzero orbital 
magnetization {J^j Tj) > 0, together with the stag- 
gered term in the ferro-orbital bond operator Oj, imply 
that nn even bonds are different than odd bonds. As a 
result, the model has translational invariance only by two 
sites. Only at /i = is the Hamiltonian invariant under 
the transformation that combines translation by one site 
and orbital inversion — > — r^. 

For J2 >• Ji, J(,A, we expect the spin chain to have 
a finite Haldane gap, close to the value for decoupled 
Haldane chains « O.4IJ2. Assuming that this is the 
case, we take the expectation value of spin operators in 
the ground state in order to write down the effective 
orbital model. The spin correlations are strongest be- 
tween nnn. Between nn, the correlation is weaker and 
oscillates due to the field-induced orbital dimerization, 
(Sj • Sj+i) = C[l + (-1)^'^], with C - J1/J2 < 1. The 
transverse orbital couplings are generated at second or- 
der in the spin-orbit coupling. The low energy effective 
orbital model becomes 

Hfz = V{A[l+C(-l)^]r/r/+i-ftr;+«;2rfrf , J, (13) 



j 



where A A + C(Ji + 2J(), C = CA{Ji + 2j;)/A < 1, 
h = h- CAJi, and K2 ~ AVA,. In Eq. im we have 
neglected orbital interactions beyond nnn as well as terms 
of higher order in the inverse spin gap. Among these is 
the nn transverse coupling, kit^t^j^-^, with ki ^ CA^/Ag. 



For K2 = 0, Eq. ([13| reduces to a classical Ising chain 
with a critical point at A = /i separating the AFOz from 
the PO state as long as \C,\ < 1. For any value of |C| < 
1, at the critical point all classical states satisfying the 
constraint of no nn down pseudospins are degenerate. 
Setting C = 0, we obtain -ffoff,2 in Eq. ([3| with coupling 



constants (A,/i, k) — >■ (A, /i, K2). Note that by setting 
Ki = in model (13 1 we have artificially (i) restored the 



translation symmetry by one site and (ii) enlarged the 
Ising symmetry from Z2 to Z2<8)Z2, leading to two Ising 
order parameters 



'2j+m/ 



(14) 



where m = 0, 1 is the index for even and odd sublattices, 
respectively. However this enlarged symmetry will not 
change the phase diagram in the region of interest near 
the MCP, namely for K2 < A w /i. 



IV. EXACT SOLUTION OF A U(l)(g)U(l) 
SYMMETRIC XXZ MODEL 

To our knowledge, the phase diagram of i/eff'.2 ^SiS not 
been studied. In order to examine the possibility of a 
critical phase we again consider the deformation of the 
model in Eq. ([5|, this time with ^ = 2, by perturbing 
around a U(1)(8)U(1) symmetric model. The point of in- 
terest is e = 1, but we start from e ^ 1. 

In terms of JW fermions, the model with e = near 
the MCP describes fermions hopping along the legs of a 
zigzag ladder. 



H 



XXZ, 2 



L 

E 



[4' 



H.c.) - (A 



/i)c]cj 



(15) 



For A/k — )■ 00, the strong nn repulsion can be replaced 
by the local constraint fijhj+i — for all j; see Fig. [2] 

Due to the constraint and the absence of nn hopping 
(ki = 0), any fermion configuration can be classified ac- 
cording to alternating domains of fermions hopping on 
even and odd sublattices: starting from the left end of 
an open chain, we find first rii fermions on the even sub- 
lattice, then n2 fermions on the odd sublattice, then 713 
fermions on the even sublattice, and so on (for definite- 
ness we assume ni 7^ for states with nonzero fermion 
filling V — N/L, and L even). Importantly, besides the 
trivial conservation of the number of fermions on even 
and odd sublattices, we have a set of good quantum num- 
bers Tia, a — 1,2, ...Nd, where Nd is the number of do- 



hopping 

infinite repulsion 



Figure 2: Schematic representation of the exactly solvable 
U(1)(8)U(1) symmetric zigzag model with A/k — >■ 00, corre- 
sponding to spinless fermions hopping on either even or odd 
sublattices, with infinite repulsion along zigzag bonds. 
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mains for open boundary conditions. We now define ficti- 
tious fermions with annihilation operator a^' that act on 
the Hilbert space of a fictitious lattice of reduced length. 
Let ra^i denote the position of the i-th (1 < i < na) 
fermion in the a-th domain in the original zigzag chain, 
measured from the left. This corresponds to the i'-th 
fermion. z' = X]a'<Q'^"' "I"*' zigzag chain. The 

reduced lattice is defined such that the position of the i- 
th fermion in the a-th domain in the original lattice gets 
mapped to 



fermions 



l)/2. 



(16) 



The length of the fictitious lattice is L' = (L - iVrf + l)/2 
for Nd odd or L' ^ {L - Nd)/2 for even. For an 
example of a configuration on the original and fictitious 
lattices, see Fig. [3] 

It is the Pauli exclusion for the a^' fermions that se- 
lects, after the prescribed mapping, only the states which 
satisfy the nn constraint in the original system. Accord- 
ingly, we find that the effective Hamiltonian in the low 
energy subspace is noninteracting: 



L' 

Hxxz,2 = 5][^(a]a^.+i+H.c.)-(A-/.)a]c 



A/k 



oo. 



(17) 



Therefore, the ground state energy can be promptly eval- 
uated as a function of fermion density v and domain wall 
density q = {Nd ~ 1)/^; we find 



-(l-,)sm— . (18) 



Minimizing with respect to v at fixed h. we find that for 
|/i — A| < k/2 the ground state has a finite density of 



a) 



b) 




Figure 3: a) Example of a configuration with Nd = 2 domains 
in the original zigzag chain, with n\ =2 fermions in the first 
domain at positions r\^\ = 2, r\^i = 6, and with n2 = 1 
fermion in the second domain located at position r2,i = 9. 
b) At the fictitious lattice with length L' = 4, constructed 
by keeping only the sites inside dashed rectangles in (a), the 
occupied sites are = \, r'2 = 3, and r'-^ = 4. In the original 
lattice there are only two possible hopping processes denoted 
by arcs without crosses, compatible with the infinite repul- 
sion along zigzag bonds. The interaction forbidden hopping 
process are blocked in the fictitious lattice due to the Pauli 
principle. 



1 /i- A 

vq = — arccos — 



(19) 



but zero density of domain walls. The chemical potential 
for domain walls 

[Id = {degs/dq)\^^ = (K/47r) sin 27ri/o ~ {h - A)vo (20) 

vanishes only for i/q — 0. Thus the c fermions prefer to oc- 
cupy a single sublattice, giving a finite value of 77^. There- 
fore, the ground state of Hxxz,2 in the limit A/k — > 00 
for — A| < K/2isa LL that breaks bond refiection sym- 
metry. Only intrachain correlations within the partially 
occupied sublattice decay algebraically. Notice that the 
nn Ki process mentioned below Eq. ( 13 ) involves interleg 
hopping, which creates domain walls. Due to the finite 
energy cost for domain walls, the liquid phase is stable 
against small ki perturbations. 

The phase diagram of Hxxz.2 near the MCP is plot- 
ted in Fig. [4|l, and the symmetry properties of the three 
phases that appear around the MCP are given in Table [iTl 
The LL phase with finite 77^ occurring for \h ~ A\ < k/2 
is denoted OLz. Its cartoon picture corresponds to pseu- 
dospins polarized along the field (f) in one sublattice and 
fiuctuating in the xy plane for the other sublattice; see 
Table |llj For h ~ A > k/2, Eq. ([l8| is minimized at 
i>o — 0, giving the PO phase with its exact classical 
ground state shown in Table |lT] For h ~ A < —k/2, 
Eq. (18 1 is minimized at i'q = 1/2, giving the AFOz 



phase. 

Next, we consider the effect of a weak perturbation 
e <C 1 in Eq. ([5| with respect to the LL fixed points in 
the OLz phase. Since domain walls have a finite energy 
cost fid, they are not created at small enough e. Thus 
{fi2j+m) — for one sublattice, m = or 1 (hence rj^ ^ Q 
still holds). Setting {T2j+„i) = 1/2 in Eq. (Isj), we obtain 
an exactly solvable anisotropic XY model lor the other 



IK/A 




b) 



♦ K/A 




hl^-\ 



hlA-\ 



Figure 4: a) Phase diagram of the U(1)(8>U(1) symmetric 
model //xxz,2 in Eq. The symmetry properties of the 

phases are given in table IIII b) Adding an infinitesimal e per- 
turbation in Hc{t,2 in Eq. Tsf leads to the opening of an energy 
gap in the OLz phase of a); the resulting phase AFOyz has 
coexisting order parameters as specified in table [III 
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2 (gapless) 




AF0y4 
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Table II: Symmetry characterization of the phases stabilized 
by the effective orbital model in Eq. (|5| with £ = 2. For 
an intuitive picture we include cartoons of the orbital con- 
figurations, with It"" = |) = I t), = -|) = I i), 
= I) = I = -|) = I and |.) denotes a 
fiuctuating orbital. The phases are characterized in terms of 
the order parameters in Eqs. Q and ( |14[ ). 



m = 1 — m sublattice, 



XY 



E 



w(l +e) y y 

2 '''2j+ni''2j+2+m 



(21) 



with e — 5- 0^. It is readily found that an energy gap 
opens for finite e in the region previously occupied by 
the gapless OLz phase. Physically, the e perturbation 
lowers the U(1)(8)U(1) symmetry down to 1i®T^2- The 
Ising Z2 symmetry acting on the fluctuating sublattice 
TO of the OLz phase is now broken, resulting in a phase 
with 7^ 0. This result also follows from the field the- 
ory analysis of the instability of the LL under the U(l) 
symmetry breaking perturbation discussed in Sec. |II A[ 
it can be verified by taking derivatives of the ground state 
energy that the Luttinger parameter of the OLz phase is 
K = \ (free fermions), independent of v. As a result, 
U(l) symmetry breaking is a relevant perturbation for 
all values of v. This is in contrast with the result for 
flxxz.ip^the crucial difference is that the size of the re- 
duced lattice defined here for the zigzag chain does not 
depend on the number of fermions, only on the number 
of domains, which is always Nd = 1 for the ground state. 
The resulting phase diagram at fixed infinitesimal e is 
shown in Fig. ^p. The phase that replaces the OLz is 
denoted AFOyz due to the coexisting 77^ and ryj^ order 
parameters; see also Table ITT] 



The effective XY model (211 applies also along the 
transition from AFOyz to the AFOz phase. Physically, 
even at finite e we expect the AFOz- AFOyz transition to 
be of Ising type and to occur within one sublattice, with 
the other sublattice being essentially frozen by the field. 
However, the transition to the PO phase, at which the 
gap for the creation of domain walls vanishes, /i^ -> 0, 
needs further consideration at finite e. For this reason, 
we shall now expand our treatment of the model to the 
global phase diagram, based on general field theory ar- 
guments. 



V. GLOBAL PHASE DIAGRAM OF Za^Za 
COUPLED ISING CHAINS 

In this section we argue that the global phase diagram 
of 7?cff,2 in Eq. ([3| has the form shown in Fig. [sj One of 
our motivations is to understand whether there is a direct 
transition between the PO and AFOyz phases near the 
MCP. 

The main argument is based on the effective field the- 
ory starting in the limit k ^ ft, 3> A. For n/h finite 
and h/A — ^ 00, -ffeff,2 in Eq. ([3| is equivalent to two 
decoupled transverse field Ising chains on the even and 
odd sublattices. At large we expect a fourfold de- 
generate phase, denoted AF0y4, which breaks both Z2 
Ising symmetries; see Table [TT] If we decrease the ratio 
k/H keeping h/A large, there is a transition from this 
ordered phase to the PO phase on a line given asymptot- 
ically by h — k/2. Point O in Fig. [s] represents the limit 
h = k/2, k/A 00. At this point the AF0y4-P0 tran- 
sition is described by two decoupled Majorana fermion 
theories with Hamiltonian densitv^S 



p=i 

(22) 

with p the sublattice index, velocity vq — k/2 and masses 
TOl — TO2 — h ~ k/2. 

However, the theory of the generic AF0y4-P0 tran- 
sition, % — Ho + Ha, includes a correction due to the 
coupling oc A between the two Ising chains. Since in the 
continuum 



we obtain 



T2X+P - iXpLix)XpRix), 



Ha = -2AxiLXiflX2LX2fl,- 



(23) 



(24) 



Defining one Dirac fermion with left and right moving 
components. 



L,R 



XlL,B. + ^X2L,B 

V2 



(25) 



we may decompose the operator Ha — y (V'l ^ 



■(/'|j)(V'L — V'l)(V'fi ~ V'fl) into a list of operators. At small 
A/k the most important contribution is the marginal at- 
tractive interaction ^ —^^^Lip^p^R. Consequently, the 
AF0y4-P0 transition becomes a LL theory with contin- 
uously varying exponents, as encountered in the Ashkin- 
Teller model. Upon increasing A/k further, the first op- 
erator that becomes relevant is ~ '4'l'4'l'4'R'4'R- Bosoniz- 
g-iv^(0±e) obtain 



ing^"^ with ^pL,R = 



1 

/2Tr ~ 



+go cos(v47r6') -I- g„ cos{Ay/T:(f 



(26) 
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with go cx mi = TO2, gu 



renormalized velocity 
V — vo + 0{A'^), and K I— [We ignored irrelevant 
terms like ^I/«^L/fl^fl,/LVfl,/L - 9,(0 ± 0)e-'^^(^Te) 



and VlV'lV'flV'fl 
mula used here differs from the one in Eq. ( 11 ) by a du- 
ality transformation (jj. The operator has scaling 
dimension l/K and is relevant for K > 1/2. In this rep- 
resentation, relevant go > describes the massive AF0y4 
phase and relevant go < the PO phase. The critical line 
approaching O is defined by setting 5o — 0. 

The (7„ operator in Eq. ( 26 1 has scaling dimension 
AK. As we depart from point O along the critical line, 
the LL parameter K decreases below unity. At the 
point where K reaches the value 1/2, denoted by A in 
Fig. [5] the perturbation becomes relevant and the 
system undergoes a Kosterlitz-Thouless transition into 
a gapped phase. To identify this gapped phase, recall 
that both Z2 Ising orders are on their instability along 
the OA line, mi = m2 = 0. Upon crossing point A 
the field (j) gets locked in the minima of the cosine 
term, hence the order parameter [using Eqs. (23), (25) 
and ^Pl,r = ^e-'v^('^±«)[ 



Notice the bosonization for- 



V'' ~ iXiLXiR - iX2LX2R ~ sin(V 47r0) (27) 

becomes finite. This implies opposite expectation values 
for the two mass terms in T-Lq. Equivalently, within a 




hl^ 



Figure 5: Schematic global phase diagram of model (|3| with 
1 = 2, namely with NNN transverse couplings. All phases are 
gapped and characterized according to the broken symmetries 
as given in Table [IT] Starting from the disordered PO phase, 
both Ising order parameters rj^ and rj^ acquire finite values 
upon crossing line OA. In the AFOyz phase the r;^ order pa- 
rameter is finite, as well as one of the two r)^ (m = 0, 1) or- 
der parameters. All three critical lines emanating from point 
A are described by a Luttinger liquid theory, which follows 
from a formal analogy between the effective theory of point 
A, Eq. ( |26[ ), and the effective theory of the XYZ spin chain 
Eq. (28 1. Point A has an emergent SU(2) symmetry, corre- 



sponding to the Heisenberg point of the XYZ spin chai n, an d 
realized in terms of the dual ^ variables defined in Eq. (29 1. 



simple mean field picture valid at finite r]^ , the interaction 
T~La gives rise to a finite mass difference mi — m2 oc 77^. 
As one mass term becomes positive and the other one 
negative, one Ising Z2 symmetry is spontaneously broken, 
leading to rf^ ^ for m — or 1. Therefore, upon 
crossing A along OA, we reencounter the AFOyz phase 
discussed in Sec. II VI 



Interestingly, the field theory in Eq. ( 26 1 controlling 



point A coincides with the field theory of the XYZ model 



H 



XYZ 



E 



(28) 



near the Heisenberg point Jx = Jy = Jz.'^In the XYZ 
model three gapped Neel phases with staggered magneti- 
zation along a = x,y ov z meet at the Heisenberg point, 
and are separated by LL transitions occurring when the 
two largest Ja's are equal (see also Ref. [50]) . From this 
field theory equivalence, we conclude that the LL transi- 
tion AO splits at point A into two LL transitions; see Fig. 
|5] It is natural to infer a LL type AFOyz-PO transition 
extending all the way to the MCP, connecting points A 
and B as in Fig. [5] 

Defining dual variables (normalized as Pauli matrices) , 



2-/ = 



2rj'. 



(29) 



in terms of which the model becomes 
1 



H, 



cff,2 



4E( 

3 



(Am^-M^\2 + A^Mj/i.^+i + 2V-M,\i), (30) 



we see that our model is self-dual along the line /i = 0. As 
we expect a single transition along the k axis, separating 
AFOz and AF0y4, this point denoted C in Fig. |5]must 
occur at K = A. It is natural to connect this point with 
the third LL line coming out of A. Finally, we expect an 
Ising transition extending from C to B and separating 
the AFOz and AFOyz phases, which differ by a single 
Ising order parameter. 

We mention that bosonizing the dual model using 
/i| ^ -^dx<p we obtain back Eq. (26). In terms of dual 
variables the AO line is an XY model, and the perpen- 
dicular direction corresponds to finite Jx — Jy At point 
A the g^ perturbation, which corresponds to the umk- 
lapp operator in spin chain languagej'^Zl becomes relevant. 
The value K — 1/2 dX point A implies an emergent SU(2) 
symmetr jffS manifested in the correlation function for the 
components of the /x vector. 

The global phase diagram changes upon including 
small finite nn transverse coupling ki which breaks the 
Ising 'Li®'Li symmetry down to Z2. Most notably the 
degeneracy of the AF0y4 phase is lowered down to 2, 
and we expect AC and AO to become Ising type transi- 
tions. On the other hand the degeneracies of the three 
phases occurring near point B are unchanged and we ex- 
pect the structure of the phase diagram in the vicinity 
of this MCP to remain unmodified up to some finite k,\. 
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Recall that in the limit K2 ^ ki ^ A the phase diagram 
is the one in Fig. [T] We leave the evolution of the phase 
diagram as a function of Ki/k2 as an interesting open 
question. 

VI. DISCUSSION 

In summary, we explored orbital ground states in the 
vicinity of a macroscopically degenerate classical transi- 
tion. We have explicitly discussed ID spin gapped sys- 
tems where the energy scale of the orbital excitations is 
pushed down to low energies due to the competition be- 
tween the Ising type interaction and the energy splitting 
induced by a lattice distortion (orbital field). Near the 
critical point, quantum fluctuations generated by spin- 
orbit coupling become important and turn the classical 
critical point into a multicritical point. 

For the toy model in which a classical Ising orbital 
chain and a Haldane spin chain are coupled by weak rel- 
ativistic spin-orbit interaction, we find that the effective 
orbital model has a multicritical point in the universality 
class of the quantum ID ANNNI model. Besides the clas- 
sical para-orbital and antiferro-orbital ordered phases, 
there appear near the multicritical point a phase with 
transverse orbital order and a critical orbital liquid phase. 
This Luttinger liquid phase arises without continuous or- 
bital symmetry in the lattice model, in the region of the 
phase diagram between ordered phases with independent 
order parameters. 

Motivated by experimental results for the compound 
CaV204, we also considered the effective orbital model 



for a zigzag chain in the regime of dominant next-nearest- 
neighbor exchange coupling. We have found an exact 
solution for a U(1)(8)U(1) symmetric model equivalent to 
hard-core fermions with hopping along the legs of the 
zigzag chain and a nearest-neighbor exclusion constraint. 
This model contains a doubly degenerate critical phase 
where the Luttinger liquid is confined to one of the legs. 

More realistically, in the model where the U(1)(8)U(1) 
symmetry is broken down to Z2 (which stems from time 
reversal symmetry in the original spin-orbital model) , we 
predict that the Luttinger liquid phase does not survive 
near the multicritical point, but is replaced by a gapped 
phase with coexisting orbital order parameters. The in- 
terpretation for this phase is that the two legs have differ- 
ent polarization in the direction of the orbital field, as in 
the classical ordered phase, but at the same time the Z2 
symmetry is broken. This result should be relevant for 
CaV204 if the appropriate model parameters satisfy the 
hierarchy of energy scales assumed here, namely that the 
next-nearest-neighbor spin exchange is the dominant en- 
ergy scale, much larger than the spin-orbit interaction. If 
the compound is also in the regime where the orbital field 
is comparable to the orbital Ising interaction, it would be 
interesting to search experimentally for the orbital mul- 
ticritical point. 
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